Modeling of Tunneling Spectroscopy in HTSC 
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The tunneling density of states of HTSC is calculated taking into account tight binding band 
structure, group velocity, and tunneling directionality for s-wave and d-wave gap symmetry. The 
characteristic density of states has quasiparticle peaks' asymmetry, flat s-wave and cusplike d-wave 
case subgap behavior, and asymmetric background. We consider that the underlying asymmetry of 
the conductance peaks is primarily due to the features of quasiparticle energy spectrum, and the d- 
wave symmetry enhances the degree of the peaks' asymmetry. Increasing of the lifetime broadening 
factor changes the degree of tunneling conductance peaks' asymmetry, and leads to the confluence 
of the quasiparticle and van Hove singularity peaks. 



I. INTRODUCTION 

Tunneling measurements on HTSC have revealed a rich 
variety of properties and characteristics [1-4] . They may 
be classified according to their low and high energy fea- 
tures. With the low energy features we may attribute: 
(i) variable subgap shape of conductance, ranging from 
sharp, cusplike, to a flat, BCS-like feature [1]; (ii) volt- 
age and temperature dependence of quasiparticle con- 
ductivity [5,6]; (iii) subgap structure [2]; (iv) zero bias 
conductance peak (ZBCP) [7]. The high energy features 
include: (i) asymmetry of conductance peaks [1] , (ii) van 
Hove singularity (VHS), (iii) conductance shape outside 
of the gap region (background (BG)) and its asymme- 
try [1]; (iv) dip feature [8]; (v) hump feature [8]. These 
features are collected in schematic Fig 1 . 
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While the tunneling spectroscopy on conventional su- 
perconductors allows directly to find the energy gap of 
superconductor, the same measurements in HTSC are 
not as easily interpreted. Some times the same exper- 
iments on the same samples show different results [9] : 
cusplike or flat subgap feature, symmetric or asymmet- 



ric conductance peaks. Usually the sharpest gap fea- 
tures are obtained when BG is weakly decreasing. A 
quantitative measure of it is the ratio of the conduc- 
tance peak height (PH) to the background conductance: 
PHB=PH/BG. When the BG conductance is decreasing, 
the PHB > 2, but when BG conductance is linearly 
increasing (~V), PHB < 2 . Kouznetsov and Coffey 
[10] and Kirtly and Scalapino [11] suggested that the lin- 
early increasing BG is arising from inelastic tunneling. 
As was mentioned in [1], the conductance is dominated 
by quasiparticle tunneling and that the effect of Andreev 
reflection is not significant. A theoretical model for tun- 
neling spectroscopy employing tight-binding band struc- 
ture, d x 2_ y 2 gap symmetry, group velocity and tunneling 
directionality was studied by Z. Yusof, J. F. Zasadzinski, 
L. Coffey and N. Miyahawa [1]. An angle resolved pho- 
toemission spectroscopy (ARPES) band structure spe- 
cific to optimally-doped BSCCO (Bi-2212) was used to 
calculate the tunneling density of states for a direct com- 
parison to the experimental tunneling conductance. This 
model produces an asymmetric, decreasing conductance 
background, asymmetric conductance peaks and variable 
subgap shape, ranging from sharp, cusplike to a flat, 
BCS-like feature. A standard technique in analyzing the 
tunneling conductance is to use a smeared BCS function 



N(E) = N(0) 



E-iT 



y/(E - iTf - A 2 



(1) 



in which a scattering rate parameter (lifetime broaden- 
ing factor) r is used to take into account any broadening 
of the gap region in the DOS. Fig. 2 shows the DOS 
calculated by formula (1) at A = 46 meV and T = 9 
meV (a), T = 3 meV (b) and T = (c). The charac- 
teristic features of the DOS is the flat subgap structure 
at small I\ This method can not explain the asymme- 
try of the conductance peaks observed in the tunneling 
experiments. 

In the case of d-wave symmetry we have 
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E-iT 
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and DOS calculated by this formula are presented in Fig. 
3. The characteristic features of the DOS is the cusplikc 
subgap structure. As was mentioned in [1] this stan- 
dard technique requires that the comparison be made 
with normalized tunneling conductance date, and since 
HTSC tunneling conductance can exhibit a varied and 
complex background shape, this procedure may "filter 
out" too much information from the conductance data. 
An alternative is to simply normalize the data by a con- 
stant. 

In [8] the tunneling data were first normalized by con- 
structing a "normal state" conductance obtained by fit- 
ting the high bias data to a third order polynomial. 
The normalized conductance date were compared to a 
weighted momentum averaged d-wave DOS 

Here /(</>) is an angular weighting function, which al- 
lows for a better fit with the experimental date in the 
gap region. A weighting function f(<f>) = 1 + 0.4cos(4</>) 
was used which imposed a preferential angular selection 
of the DOS along the absolute maximum of the d-wave 
gap and tapers off towards the nodes of the gap. This 
is a rather weak directional function since the minimum 
of /(0) along the nodes of the d-wave gap is still none- 
negligible [8]. 

A. J. Fedro and D. Koclling [12] have done the mod- 
eling of the normal state and superconducting DOS of 
HTSC, using tight-binding band structure, including the 
next nearest neighbors 



£fc = —2t(cos(k x a) + cos(k y aj) + 4t' cos(k x a)cos(k v a) — /i 

(4) 



The calculation showed two singularities in DOS: a 
van Hove singularity in the center of energy band due 
to saddle point near (n,0) at t' = and another at the 
lower edge of the energy band due to extra flattening out 
at (0,0). As extended s-wave and d-wave superconduct- 
ing DOS were considered in case of hole-doped situation 
(/j < 0) for different hole concentration. The Fermi sur- 
face for t' = and t' = 0.45t at the same concentration, 
corresponding [i/2t = —0.187 which was used in [12] are 
presented in Fig. 4a. It is needed to move up from the 
Fermi surface (set as the zero energy) to reach the point 
(7r,0) in case t' = and move down in case t' = 0.45t. 
So, for t' — the Fermi energy lies to the left of the 
van Hove singularity and will move away from it with 
increased hole doping while for t' = 0.45i it lies to the 
right and will move towards to with increased hole dop- 
ing (See Fig. 4b where the DOS for t' = and t' = 0.45 
are presented ) . For calculation of the superconducting 
DOS Fedro and Koelling used formula 

N(E) = \ £(1 + ^-)S(E E k ) + (1 - ^-)6(E + E k ) 



Ek ' 



(5) 



This formula is the limit of the expression for tunnel- 
ing density of states (6) at T — and \T k \ 2 = 1, where 
Tfe is tunneling matrix element. 
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Fig. 4 Fermi surfaces (left) and DOS (right) for t'—O (solid lines) and i'=0.45t (dashed lines) in formula (4) at 



——-0.187 which corresponds hole - doped situation. 
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Fig. 5 DOS for t'=0 at different T for s-wave symmetry (a) and d-wave symmetry (b), calculated by formula (6). 



The Fig. 5a shows the result of calculation of the DOS 
for t' = at r = 0.07, 0.1 and 0.2 meV for s-wave sym- 
metry which reflect the results of Fedro and Koclling. 
Fig. 5b shows the same DOS for d-wave symmetry. In 
both cases the Fermi energy (set as the zero of energy) 
lies to the left of the van Hove singularity. There is the 
peaks' asymmetry which is more pronounced at large T. 



ul = - 2 (l 
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V * = 2 (1 - E-J 



(8) 



and r is the quasiparticle lifetime broadening factor. The 
energy spectrum of quasiparticles in the superconducting 
state is determined by 



II. MODELS AND METHODS 

In this paper we use the method for calculation of the 
DOS presented in [1]. The tunneling DOS of a super- 
conductor is determined by the imaginary part of the 
retarded single particle Green's function 



N(E) = --I m y2\T k \ 2 G R (k,E) 

k 

For the superconducting state 

„.2 „,2 

G R (k,E) 



E-E k +iT E + E k 



where uj. and v\ are the usual coherence factors, 



(6) 



(7) 



E, 



+e k 



(9) 



with the effective band structure extracted from ARPES 
experiments [13] 

6c = Co + 0.5Ci[cos(k x a) + cos(k y a)] 
+ C'2Cos(k x a)cos(k y a) + Q.'bCz[cos{2k x a) + cos(2k y a)] 
+ 0.5C4[cos(2k x a)cos(k y a) + cos(k x a)cos(2k y a)] 

+C^cos(2k x a)cos(2k y a) (10) 

Here ^fc is measured with respect to the Fermi energy 
(6=0), and the phenomenological parameters are (in 
units of eV) C = 0.1305, d = -0.5951, C 2 = 0.1636, 
C* 3 = -0.0519, C 4 = -0.1117, C 5 = 0.0510. 
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Fig.6(left) 3D-plot of energy spectrum of normal state according to formula (9). 
Fig.7(center) 3D-plot of coherence factor u\ according to formula (8). 

Fig.8(right) Fermi surface corresponding to the = in formula (10). Dark straight line shows the line of directional 
tunneling, the dashed lines show the angular spread Oq. 



Fig. 6 shows the three dimensional image of function 
(10) . There are saddle point in (nfl) and flattening out 
of the energy band at (0,0) which lead to the van Hove 
singularities in the DOS. The three dimensional graph of 
the coherence factor u\ is shown in Fig. 7. 

Since quasiparticles with momentum perpendicular to 
the barrier interface have the highest probability of tun- 
neling, the tunneling matrix element \Tk\ 2 reveals a need 
for factors of directionality D(k) and group velocity v g (k) 
[1]. The group velocity factor is defined by 

v g (k) = |V fc &.n| = \§^cos(6) + §rsin(6)\ (11) 

where the unit vector n defines the tunneling direction 
as shown in Fig. 8, which is perpendicular to the plane 
of the junction. 

The directionality function D(k) is defined by 

Here 0o defines the angular spread of the quasiparti- 
cle momentum with none-negligible tunneling probabil- 
ity with respect to n. The tunneling matrix element |7fc| 2 
is written as 



\T k \ 2 = v g (k)D(k) (13) 

The three dimensional graphs of group velocity v g (k), 
directionality D(k) and tunneling matrix element \Tk\ 2 
functions are shown in Fig 9. 



III. RESULTS AND DISCUSSIONS 



Different factors may lead to the changing of the en- 
ergy gap A in HTSC. In particular, strong effects are 
caused by nonmagnetic impurities [14]. In superconduc- 
tors with d-wave symmetry the nonmagnetic impurities 
destroy the superconductivity very efficiently. Possibil- 
ity to destroy of Cooper pairs by impurities leads to their 
finite lifetime. If the state with the quasiparticle is not 
stationary state, it must attenuate with time due to tran- 
sitions to other states. The corresponding wave function 
has the form e~ 4 ^ p )*/ fi ~ r */ ?i , where V is proportional to 
the probability of the transitions to the other states. It 
may be interpreted as the energy of the quasiparticle has 
the imaginary addition — iT. The relation between Y and 
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Fig. 10 The changing of DOS with energy gap Ao for s- (a,b) and d-wave (c,d) symmetry at T—3 meV (a,c) and 
T=9 mcV (b,d) without of effects of directionality and group velocity. 



lifetime of quasiparticlc r s is T = h/r s . Hence, the im- 
purities lead to changing Ao and we may do modeling of 
the influence of impurities on tunneling conductance by 
numerical calculation of DOS N(E) considering different 
values of Ao in formula (6). Here we present results of 
calculation N(E) at A = aA o, where a— 0.2, 0.4, 0.6, 
0.8, 1 and A 00 = 46 mcV. 

The peculiarities of the quasiparticle energy spectrum 
(10) play an essential role in explanation of the conduc- 
tance features. Here, based on the numerical calculation 
of DOS we consider that the underlying asymmetry of 
the conductance peaks is primarily due to the features of 
quasiparticles energy spectrum. The d-wave gap symme- 
try simply enhances the degree of the peaks asymmetry. 
The last one is also changed by changing the tunneling 
direction. 

Fig. 10 shows the results of the numerical calculations 
of the DOS at T = 3 mcV (a,c) and r = 9 mcV (b,d) 
without effects of group velocity and directionality as for 
s-wave (a,b) and d-wave (c,d) gap symmetry, correspond- 
ingly for different values of energy gap Ao- We have de- 
creased the energy gap A , starting from A = 46meV. 
For clarity we present only three characteristic curves, 
which corresponds to aAo with a =1, 0.6 and 0.2. We 
exclude the effects of group velocity and directionality 
to demonstrate that they are not responsible for peaks 
asymmetry. 



There is the asymmetry of the quasiparticle peak 
heights as for s- and d-wave symmetry. So, the origin 
of the peaks asymmetry is not due to d-wave symmetry 
of the energy gap of HTSC. There is more flat subgap 
behavior of DOS in the case of s-wave symmetry in com- 
paring with the d-wave case. The increase of lifetime 
broadening factor T leads to the enhance of the peaks' 
asymmetry. There are van Hove singularities in the DOS 
at small T. The increase of T leads to the confluence 
of the quasiparticle and VHS peaks and this results to 
the enhance of the DOS peaks asymmetry due to saddle 
point in energy spectrum (10) at (7T,0). Also note to the 
asymmetry of the background as for s- and d-wave gap 
symmetry. 

Fig. 11 shows the Ao-dependence of DOS taking into 
account the effects of group velocity and directionality at 
T = 3 meV (a,c) and T = 9 mcV (b,d) for s-wave (a,b) 
and d-wave (c,d) gap symmetry. As in [1] we have taken 
= 0.25 and 6 = 0.1. 

There is also quasiparticle peaks' asymmetry similar to 
s- and d-wave cases. But in d-wave case the asymmetry 
is more stronger than in s-wave. The effects of group ve- 
locity and directionality lead to disappear of the VHS in 
DOS. The increase of T enhances the quasiparticle peak 
asymmetry. The most strong effect of energy band struc- 
ture on the DOS occurs along A: x -axis due to van Hove 
singularity at (tt,0). 
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Fig. 11 The changing of DOS with energy gap A for s- wave (a,b) and d-wave (c,d) gap symmetry at T=3 meV 
(a,c) and T=9 meV (b,d) with the effects of directionality and group velocity. 
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Fig. 12 Effects of directionality on the DOS as in s-wave (a,b) and d-wave gap symmetry at T—3 meV (a,c) and 
r=9 meV (b,d). 
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Fig. 12 has demonstrated this effect. We have pre- 
sented the DOS at different © at T — 3 meV (a,c) and 
T = 9 meV (b,d) as for s-wave (a,b) and d-wave (c,d) 
gap symmetry. In the case of s-symmetry the position of 
the quasiparticle peaks is constant excluding the direc- 
tion along k x { 6 = 0). We pay attention to the strong 
peaks' asymmetry in this case. 

In the case of d-wave symmetry we have practically 
the same behavior around k x direction as for s-wave, but 
energy gap is changed due to the 0-dependence of A 
and correspondingly, the quasiparticle peaks are shifted 
to the zero energy. 

Fig. 13 shows the 6o-changing of DOS at Y = 3 meV 
(a,c) and r = 9 meV (b,d) as for s-wave (a,b) and d-wave 
(c,d) gap symmetry. The increase of ©o means the tak- 
ing into play (inclusion) the states, close to (7r,0). It is 
reflected as an appearance of the van Hove singularity as 
in case of s-wave and d-wave gap symmetry at small V. 
The VHS is more pronounced in case of d-wave in com- 
paring with s-wave symmetry. The increase of Y leads to 
confluence of the quasiparticle and VHS peaks. 

We consider that the absence VHS peak on the ex- 
perimental <AZy<iy-characteristics means the enough large 
lifetime broadening factor V in that HTSC material. 

The origin of the peaks asymmetry in the tunneling 
DOS was studied in [1] by considering the role of the 



tunneling matrix element \Tk\ 2 in the clean limit case 
(r = 0), where for the calculation N(E) was used the 
formula (5). 

We repeat the explanation of the paper [1] because we 
believe that the following conclusion on the origin of the 
peaks asymmetry must be different. At E > (positive 
bias voltages) the first term of (5) contributes to the N(E) 
because of 8(Ek — E). In this case, as can see from Fig. 
7 and Fig. 9c and Fig. 14a \Tk\ 2 selects only a relatively 
short region of states in k-space in which uj, > 0. These 
are the states with > (above the FS). For the major- 
ity of states integrated over ( see again Fig. 7 and Fig. 9c). 
At E < (negative bias voltages) the second term of (5) 
contributes to the DOS because of 5(Ek + E). In this 
case, as can see from Fig. 7 and Fig. 9c and Fig. 14b, 
\Tk\ 2 selects out a large region of k states where v\ > 0, 
in fact, equal to one. These states are below the Fermi 
surface, where < 0. The overall effect then is to have a 
large negative bias conductance compared to the positive 
one. This is true as for s- and d-wave symmetry. Hence, 
the underlying asymmetry of the conductance peaks is 
primarily due to the band structure and d-wave sym- 
metry simply enhances the degree of asymmetry of the 
peaks. So, the peaks' asymmetry existing as for s-wave 
and d-wave symmetry is sensitive to band structure £fc. 
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Fig. 14 ARPES energy spectrum along 9 = 0.25. The 
E < (b). 

In summary, by changing of the energy gap Ao in 
HTSC one may model the influence of nonmagnetic im- 
purities on the DOS. We consider that the asymmetry of 
the quasiparticle peaks is due to the specific features of 
the energy spectrum of HTSC and that the d-wave gap 
symmetry only enhances the peaks' asymmetry. 

The absence of the VHS peak on the experimental 
dl/eiP-charactristics means the large enough lifetime 




values of coherent factors correspond to E > (a) and 
broadening factor T in HTSC. 
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